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fracted by the far limb, and then fitted to the data. Values of 
the pressure, derived from separate immersion and emersion 
chords, show no significant trends with latitude, indicating 
that Triton’s atmosphere is spherically symmetric at ~50-km 
altitude to within the error of the measurements; however, 
asymmetry observed in the central flash indicates the atmo- 
sphere is not homogeneous at the lowest levels probed (~20.- 
km altitude). From the average of the 1995 occultation data, 
the equivalent-isothermal temperature of the atmosphere is 
47 ± 1 K and the atmospheric pressure at 1400-km radius 
(~50-km altitude) is 1.4 ± 0.1 /d>ar. Both of these are not 
consistent with a model based on Voyager UVS and RSS obser- 
vations in 1989 (D. F. Strobel, X. Zhu, M. E. Summers, and 
M. H. Stevens, 19%, Icarus 120, 266-289). The atmospheric 
temperature from the occultation is 5 K colder than that pre- 
dicted by the model and the observed pressure is a factor of 
1.8 greater than the model. In our opinion, the disagreement 
in temperature and pressure is probably due to modeling prob- 
lems at the microbar level, since measurements at this level 
have not previously been made. Alternatively, the difference 
could be due to seasonal change in Triton’s atmospheric 

structure. © IW Academic Press 


I. INTRODUCTION 

Our knowledge of Triton’s lower atmosphere is based on 
relatively few observations. Spectra observed from Earth 
reveal the surface ices N 2 , CH 4 , CO, and CO 2 (Cruikshank 
and Silvaggio 1979, Cruikshank et al. 1984, 1993). For a 
wide range of temperatures, these species would coexist 
as vapor and ice. N 2 is the major atmospheric constituent, 
because it has the greatest vapor pressure at the tempera- 
ture of 38 ± 1 K inferred from the band shape of N 2 (Tryka 
etal. 1993) and the Voyager IRIS measurements (Conrath 
et al. 1989). Two of the species, N 2 and CH 4 , were directly 
detected as gases by the Voyager UVS solar and stellar 
occultations, while an upper limit of 0.01 was placed on 
the mixing ratio of CO (Broadfoot et al. 1989). These UVS 
occultation data also determined the N 2 temperature at an 
altitude of 475 to 675 km to be 102 ± 3 K (Krasnopolsky 
et al. 1993). Voyager 2 images show dark plumes rising 
from the surface and leveling off at an altitude of 8 km 
(Soderblom et al. 1990). Clouds and haze were detected 
below 30-km altitude (Smith et al. 1989; Pollack et al. 1990). 

The most direct probes of the lower atmospheric struc- 
ture are the two occultations of the Voyager spacecraft 
radio signal by Triton as seen from Earth (Tyler et al. 1989, 
Gurrola 1995). The phase shift in the radio signal directly 
yields the refractivity of the atmosphere as a function of 
radius, which can be converted to a number density profile 
if one knows the mix of gases. Finally, one can obtain 
temperature and pressure profiles if one assumes the atmo- 
sphere is in hydrostatic equilibrium. Since the phase shift 
of the radio signal was small in the neutral atmosphere. 


the analysis of the RSS results was based on model fitting, 
rather than numerical inversion, and the surface pressure 
was found to be 16 ± 3 p.bar (Tyler et al. 1989), later 
refined to 14 ± 1 /xbar (Gurrola 1995). Another relevant 
result of the RSS occultations is a determination of Triton’s 
surface radius, 1355 ± 7 km (Tyler et al. 1989), which can 
be compared with the value of 1352.6 ± 2.4 km derived 
from measurements of the Voyager images (Davies et al. 
1991). 

Models of the pressure and temperature as a function 
of altitude for Triton’s atmosphere make the implicit as- 
sumption that the atmosphere is substantial enough that 
the surface ice is maintained at a constant temperature 
through buffering by sublimation and condensation (Traf- 
ton 1984, Ingersoll 1990, Yelle et al. 1995). Given this 
assumption, one can use atmospheric measurements from 
different locations on Triton to constrain a single atmo- 
spheric model of the vertical temperature and pressure 
profiles. Models by Strobel and Summers (1995), Krasno- 
polsky etal. (1993), and Strobel et al. (1996) determine the 
vertical atmospheric structure above the tropopause from 
radiative-conductive calculations, which include the ef- 
fects of CH 4 heating and cooling, CO cooling, magneto- 
spheric electron heating, and solar EUV heating. Broadly 
speaking, the resulting thermal profile begins at the surface 
temperature of 38 K, drops adiabatically to 37 K at an 
altitude of 8 km (following the troposphere model of Yelle 
et al. 1991), then rises smoothly to 102 K at 600 km. 

Given our current state of knowledge of Triton’s atmo- 
sphere, we set out to (1) probe Triton’s atmosphere in the 
microbar pressure region where previous observations do 
not exist and compare with model predictions, (2) investi- 
gate seasonal change which has been predicted based on 
the vapor-pressure equilibrium of Triton’s atmosphere 
with the surface ices, and (3) investigate spatial variability 
with a multi-chord occultation observation. These tasks 
require high-resolution information about Triton s atmo- 
spheric structure, and at present this information can be 
acquired from Earth only through observation of a series 
of stellar occultations. We began such observations in 1993, 
with the successful Kuiper Airborne Observatory (KAO) 
observation of the occultation of the star Tr60, R = 13.9 
(McDonald and Elliot 1992), by Triton. On 1995 August 
14, nine additional chords were obtained from six sites 
during the occultations of Trl48A and Trl48B including 
combined visible and IR observations from two sites, ob- 
servations of both occultations at two sites, and the obser- 
vation of the central flash for Trl48B from the IRTF. These 
data probe an altitude range near the microbar pressure 
level (20-100 km of altitude) intermediate to those sensed 
by Voyager’s instruments. 

In this paper the entire ground-based stellar occultation 
data set (Tr60, Trl48A, and Trl48B) is used to address 
the goals described in the previous paragraph by fitting 
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the lightcurve model of an isothermal atmosphere (with 
an optional power-law thermal gradient) developed for 
Pluto by Elliot and Young (1992). 

II. PREDICTION OF THE GROUND TRACKS 

The occulted stars, Tr60 and Trl48, were identified in 
an observational program to search for Triton occultation 
candidates (McDonald and Elliot 1992, 1995). Since the 
error in these initial predictions was estimated to be ±0.25 
arcsec (equivalent to nearly the diameter of the Earth 
subtended at the Earth-Triton distance), further observa- 
tions of Triton and the stars were needed to improve the 
predictions, to deploy observers within the ground tracks 
of the occultations. Triton’s radius subtends an angle of 
about 0.07 arcsec as seen from the Earth, so our final 
accuracy goal, to place the KAO within the shadow with 
confidence, was ±0.02 arcsec. 

Refinement observations for occultation predictions are 
most effective when both bodies involved in the occultation 
can be recorded on the same CCD strip scan or stare frame, 
such that (unknown) differences between the coordinate 
systems of the stellar reference network and the body 
ephemeris will not affect the relative positions. Our reduc- 
tion procedure, when both bodies are on the same frame, 
allows for the transformation between the local stellar and 
ephemeris coordinate systems to include an offset, rota- 
tion, and scale factor. For a discussion of our methods used 
to reduce and analyze the prediction-refinement observa- 
tions, see Olkin (1996). 

Tr60 Prediction 

We planned three data sets that would record both bod- 
ies on the same scan or frame: (1) strip scans with the 
SNAPSHOT clone on the 0.9-m Crossley Telescope at 
Lick Observatory (Dunham etal 1985, Dunham 1995); (2) 
strip scans with SNAPSHOT (Dunham et al. 1985) on the 
0.61-m telescope at the George R. Wallace, Jr., Astrophysi- 
cal Observatory, and (3) stare frames recorded with the 
CCD camera on the 1.6-m astrometric reflector at the U.S. 
Naval Observatory (USNO) Flagstaff Station. Limited by 
the amount of sky that could be exposed on a single strip 
scan or stare frame, the Lick and Wallace observations 
began about 8 weeks before the event. The USNO observa- 
tions began just 3 days prior to the event and were used 
to make the final adjustment in the KAO flight plan. The 
final prediction based on the USNO observations gave a 
geocentric impact parameter of 0.04 ± 0.02 arcsec, which 
is within one standard deviation of the post-event value: 
0.026 ± 0.006 arcsec (see Fig. 1). 

Trl48 Prediction 

Only from the occultation observations did we learn that 
Trl48 had two components, since the two stars, Trl48A 


and Trl48B, appear blended in Earth-based imaging. Thus, 
our. prediction was based on center-of-light observations 
of the pair. With the initially predicted closest approach 
distance of 0.42 arcsec (McDonald and Elliot 1995), the 
occultation would not be visible from Earth; however, 7 
months before the occultation, astrometric observations of 
Trl48 (the blended star image) by the Carlsberg Automatic 
Meridian Circle (CAMC) yielded a geocentric closest ap- 
proach distance between the star and Triton of 0.30 ± 0.13 
arcsec (CAMC star position: J2000 a = 19 h 41 m 00?376; 8 = 
-20°51'41'/20, B. Argyle, personal communication). This 
prediction placed the centerline of the occultation shadow 
tangent to Earth’s limb with an error in the prediction 
equal to the diameter of Triton as seen from Earth. 

The prediction was further refined with both stare frame 
and strip scan observations. From 52 stare frames recorded 
at the USNO on 2 nights before the occultation, the pre- 
dicted geocentric closest approach of the center of Triton 
to the blended stellar image was 0.31 ± 0.05 arcsec. From 
strip scans recorded at Lick Observatory’s Crossley tele- 
scope (150 observations of Triton and 154 of the blended 
Trl48A-Trl48B image) the predicted closest approach 
was 0.28 ± 0.02 arcsec. For comparison, the reconstructed 
shadow path, from the multiple chord occultations (see 
Section VI), yielded geocentric impact parameters of 
0.2788 ± 0.0003 arcsec for Trl48A and 0.1961 ± 0.0012 
arcsec for Trl48B. The corresponding shadow paths are 
displayed in Fig. 2. 

To compare the predicted and reconstructed impact pa- 
rameters, we converted the reconstructed results into a 
center-of-light prediction for the pair. Using the Trl48A/ 
Trl48B intensity ratio of 4.19 (Section VIII) from unfil- 
tered visual wavelength observations of both occultations 
from one station, we calculated the center-of-light recon- 
structed prediction to be 0.247 ± 0.002 arcsec, about 1.5 
standard deviations smaller than the weighted prediction 
from the Lick and USNO data. 

III. TR60 OBSERVATIONS 

Observations of the Tr60 occultation were attempted 
from western Australia at Perth Observatory, South 
America at Cordoba Astronomical Observatory, South Af- 
rica at the South African Astronomical Observatory, and 
the South Atlantic from the KAO. A telescope at each of 
these sites was equipped with a portable CCD (PCCD) 
occultation photometer (Buie et al 1993), the first joint 
deployment of these four instruments. These observations 
were unfiltered to maximize the signal-to-noise ratio 
(SNR) of the lightcurves. Table I gives specifics of all Triton 
occultation observations recorded (or attempted) to date. 
Subframe size and detector scale are reported in binned 
pixels for those observations recorded with on-chip bin- 
ning. For datasets with lightcurves derived from aperture 
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FIG. 1. Mercator projection of Earth showing the observing station locations and a posteriori estimate of the Tr60 occultation shadow path 
which was determined from imaging data taken before and after the event at the 1.6-m USNO Flagstaff Station (geocentric impact parameter of 
0 026 ± 0.006 arcsec). The three lines of the shadow path are the northern limit, midline, and southern limit assuming a shadow radius of 1448 km. 
Australia was south of the occultation shadow. The path crossed South Africa and Argentina, but clouds prevented the observations. The event 
occurred near opposition when most of the Triton-facing hemisphere of Earth was in darkness. 


photometry, the aperture size is given; for those reduced 
using numerical point-spread function (PSF) photometry, 
the full width at half-maximum (FWHM) of the PSF is 
given. The SNR in a 1-sec integration is derived from the 
scatter in the isothermal model fit residuals. 

Since the skies in South Africa and South America were 
cloudy and the occultation path passed north of observers 
in Australia, the only observed occultation chord was re- 
corded from the Kuiper Airborne Observatory. The KAO 
was deployed from Punta Arenas, Chile, due to the ex- 
pected southern track of the occultation shadow; interme- 
diate predictions had put the shadow path as far south as 
the southern tip of South America. As we received new 
astrometric measurements, the predicted shadow path was 


refined, and each day the updated prediction consistently 
moved the shadow north. This complicated the deployment 
plans as we needed the full flight range of the KAO to get 
to the occultation shadow. 

Sample frames of the KAO observations are displayed 
in Fig. 3. We could not use aperture photometry to derive 
the lightcurve from the images because of the proximity 
of Neptune to the blended Triton-Tr60 and the rapidly 
changing image shape (due to the telescope tracking). To 
get the highest SNR lightcurve, we used a numerical PSF 
model fitting method to determine the intensity of the 
blended Triton-Tr60 image relative to Neptune. A field 
star brighter than the blended Triton-Tr60 image would 
have been the best choice as a reference PSF. Unfortu- 
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FIG. 2. Mercator projection of Earth showing the observing station locations and a posteriori estimate of the shadow path for the Trl48A 
(solid line) and Trl48B (dashed line) occultations. The occultation of both stars was recorded at two sites from Lick Observatory. Observers on 
the KAO and at WIRO recorded the occultation of the brighter star only, while observers at the IRTF and Lowell Observatory recorded the 
occultation of the fainter star only. Observations of the Trl48B occultation at WIRO and Steward Observatory were compromised by clouds, and 
those at the KAO, by instrumentation problems. The Hubble Space Telescope was south of both shadow paths. 


nately, no such star was available, so Neptune was used. 
Although Neptune is not a point source (its disk subtends 
2.3 arcsec), we can use it as a PSF source because of the 
3- to 5-arcsec seeing on the KAO. Other PSF fitting tech- 
niques, such as using an analytic model, i.e., an elliptical 
Lorentzian (Bosh et al. 1992), proved unsatisfactory be- 
cause the frequently elongated PSF was not well fit by the 
analytic model. Our method provided an error in the ratio 
for each subframe, which was used to weight the data in 
the fits to the occultation lightcurves. 

A model numerical PSF for each frame was derived 
from a cubic interpolation of the Neptune image with the 
background removed. After each iteration, the background 
was recalculated and the numerical PSF model was up- 
dated (by removing the new background). A two-source 


image model was fit simultaneously to the blended image 
of Triton-Tr60 and the image of Neptune. The fitted pa- 
rameters include an offset of the Triton-Tr60 center from 
Neptune, the ratio of the peak of the Triton-Tr60 PSF 
relative to the Neptune PSF peak, and three background 
parameters (a mean level, row slope, and column slope). 
Pixels common to both PSFs were not included in the fit. 

IV. TR148 OBSERVATIONS 

Nine lightcurves were derived from the two occultations 
in 1995. The four PCCD instruments were set up at the 
Infrared Telescope Facility on Mauna Kea (MIT system), 
aboard the KAO (NASA Ames system), at the 1.8-m tele- 
scope at Lowell Observatory (Lowell system), and at the 
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FIG. 3. Sample images from the Triton occultation of Tr60. In both images there ire two sources. The brighter object is Neptune' the other 
(lower and to the left) is the blended image of Triton-Tr60. (a) shows a “round" point-spread function (frame 1803), while (b) taken just 1 0 sec 
friT'U a Point-spread function (frame 1805). A charge transfer efficiency problem is seen by the streaking of the images toward the 

^ f FOWS ° subframe is an artifac * of ‘he detecto, readout. The image scale is approximately 1 arcsec/p.xel. 

At this time, Tnton was near eastern elongation and '"'15 arcsec from Neptune. 


Wyoming Infrared Observatory (University of Arizona 
system). At the 0.9-m Crossley telescope at Lick Observa- 
tory, a SNAPSHOT clone (Dunham et al. 1985) was used 
to record high-speed images of the occultation. This instru- 
ment, like the PCCD, was designed for high-speed occulta- 
tion observations, and the integrations are triggered from 
a GPS signal. All the visible wavelength observations were 
recorded unfiltered. The lightcurve of the blended Triton- 
Trl48 from Lowell Observatory was divided by the 
lightcurve from Neptune to remove any varying extinction. 

Two infrared instruments were also used: (1) NSFCAM 
(Shure et al. 1994) at the IRTF, and (2) LIRC2 at the 1.0- 
m telescope at Lick Observatory. The NSFCAM observa- 
tions were recorded with a K filter (center wavelength = 
2.21 /i.m, FWHM = 0.39 yum) and the LIRC2 with a K' 
filter (center wavelength = 2.12 gm, FWHM = 0.34 yum). 
For comparison, the effective wavelength of the IRTF- 
visible wavelength data is ~0.7 yum. NSFCAM has sub- 
frame readout capabilities to reduce deadtime, but does 
not have a frame transfer scheme. To establish the dead- 
time for NSFCAM, which depends on the number of pixels 
being read out and the location of those pixels, images of 
an infrared light-emitting diode blinking once per second, 
triggered by a GPS receiver, were recorded in occultation 
mode (same box size, box location, integration time, and 
series length as the occultation observations). The dead- 
time for these observations is 24 msec. With the LIRC2, full 
frame images were recorded, and the deadtime between 


integrations, while the detector was being read out, is sev- 
eral seconds. 

At WIRO, clouds of varying optical depth were present 
throughout the data recording interval. To first approxima- 
tion, the effect of clouds was removed from the blended 
Triton-Tr j 48A-Trl48B signal by using Neptune as a mon- 
itor; however, there is some evidence that the clouds are 
not completely neutral as a function of wavelength, and 
color differences between the various signal components 
may compiomise photometric accuracy at some level. 

For the KAO observations, the start time of the data 
recording was triggered from the GPS receiver, but the 
exposure timing was not, because of a communication 
problem between the instrument computer and the GPS 
receiver. As a result, the integration time was controlled 
by the clock in the instrument computer. The integration 
time was determined to be 0.342 sec from a series of images 
of the secondary mirror nodding in and out of the beam 
at a frequency of 1 Hz controlled by a GPS receiver. 

The IRTF lightcurves and their fit results differ from 
those of Olhin (1996) due to an improved method of aper- 
ture photometry used in this work. 

V. LOCATION OF THE TR60 OCCULTATION 
GROUND TRACK 

The closest approach distance between Triton and Tr60 
was best dt termined by analysis of pre- and post-event 
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imaging of the occulted star and Triton recorded at the 
1.6-m telescope at the USNO (Olkin et al 1996). These 
USNO data give the closest approach of the KAO to the 
center of Triton’s shadow to be 359 ± 133 km. As a test 
of the astrometric method, Olkin et al. (1996) found the 
occultation midtime derived from the USNO frames is 
2.2 ± 4.1 sec earlier than the midtime from the occultation 
lightcurve itself. The latitudes probed by the occultation 
lightcurve reported in Olkin et al (1996) are incorrect. 
Figure 4 depicts Triton’s globe and the path of the star as 
seen from the KAO, as well as the latitudes probed as a 
function of time. 

VI. LOCATION OF THE TR148 OCCULTATION 
GROUND TRACK 

The multiple occultation chords from the Trl48A and 
Trl48B events were used to reconstruct the occultation 
shadow path, to determine where on Triton each chord 
probed, and to establish the separation of the double star. 
First, we defined the shadow plane, also known as the fg 
plane (Elliot et al 1993), to be perpendicular to the line- 
of-sight to the star, such that the perpendicular to the fg 
plane points to the occulted star (using coordinates of the 
blended star as determined by the CAMC), /points in the 
direction of increasing right ascension, and g completes 
the right-handed system (increasing declination). An offset 
of each star from the blended center, (/ 0A , £oa) ancI (/ob* 
g 0B ), i s determined by fitting the positions of the observers 
relative to the shadow center (/ pr and g pr ) at the immersion 
and emersion half-light times to a model for Triton’s figure. 
The offsets account for errors in the star position and/or 
Triton’s ephemeris. 

The positions of the observers relative to the shadow 
center are determined by Triton’s ephemeris and the coor- 
dinates of the observers. The DE211 ephemeris and the 
NEP016 model (Jacobson et al 1991) were used to calcu- 
late a light time-corrected J2000 Triton ephemeris. The 
ephemeris was converted from right ascension, declination, 
and distance to the shadow plane coordinates / p (f) and 
g p -(0. Next, the position of the observer (receiver) in the 
shadow plane f r (t) and gft) was calculated from the geo- 
detic coordinates of the observing site (Table II). For the 
occultation of Trl48A, the position of the observer relative 
to the shadow center at a time t is given by 

/pr(0=/r(0-(/p(0-/o A ), 

£pr(0 — £r(0 ~ ( gp{0 ~ 8 oa)-> 

an analogous equation exists for the Trl48B occultation 
data. These coordinates at the immersion and emersion 
half-light times were used in a fit of Triton’s figure. The 


immersion and emersion times at half-light were derived 
from fits to the occultation lightcurves using an isothermal 
model (see Table II). The lightcurve model is an extension 
of Elliot and Young (1992) to include the flux contribution 
from the far limb and is presented in Section VII. 

The parameters of the least-squares fit include the star 
offsets and parameters describing Triton’s figure (i.e., semi- 
major axis, ellipticity, and orientation angle). These param- 
eters are adjusted to minimize the distance between the 
location probed at half-light and the figure of Triton along 
the direction of the shadow path. The orientation of the 
semimajor axis of Triton’s figure is measured from the / 
axis (positive angle measured from g to /). Because the 
lightcurves are of different quality, the fits were weighted 
with the formal errors from the half-light timings. The 
fitted parameters and their formal errors are given in Table 
III for two different models: circular figure and and ellip- 
tical figure. 

We refit the lightcurve timings to a circular figure of 
Triton without the IRTF observations (since they control 
the weighted fit) to see how robust the solution is. The 
two solutions (circular figure with and without the IRTF 
timings) agreed within their errors, and the impact parame- 
ter for the IRTF differed from the adopted solution by 
only 6 km (much less than lcr). Unlike the circular-figure 
fits, varying the data in the oblate-figure fits affects the 
results. We chose the circular-figure model as our adopted 
solution because of its robustness. The declination offsets 
(g 0A and g 0B ) of the two stars in the adopted solution are 
systematically north of the CAMC position (used as the 
reference in this system). This could be due to (1) an 
error in the Triton ephemeris, (2) an error in the CAMC 
measurement, (3) a disagreement between the coordinate 
systems of the star’s coordinate and the ephemeris, or (4) 
all the above. 

From the Faint Object Spectrometer (FOS) on the Hub- 
ble Space Telescope, we have a resolved “image” of the 
double star, exposed on the day of the occultation. Al- 
though the FOS is a one-dimensional array detector, “im- 
ages” are constructed by moving the stellar source perpen- 
dicular to the detector array to add the second dimension. 
We have fit a two-source PSF model to the image to deter- 
mine the separation of Trl48A and Trl48B. The separation 
is either greater or less than that derived from the occulta- 
tion timings (Table III), depending on the choice of refer- 
ence PSF. Therefore, these data are consistent with the 
results from the occultation timings, but cannot provide a 
more precise estimate of the star’s separation. 

From the fit results of Table III and the locations of 
the observers, we determine the impact parameter (p min ), 
which is the closest distance between the observer and the 
center of the occultation shadow (see Table IV). Also 
included in the table is the residual for each site (the dis- 
tance along the shadow path between the position probed 
at half-light and the figure) and the velocity of the occulta- 
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RG. 4. Locations on Triton probed by the stellar occultation of Tr60. (a) Path of the star at 5-sec intervals as it is occulted by Triton The 
scale bar indicates 500 km or 24 milharcsec. The uncertainty in the star path is about 6 milliarcsec. (b) Latitudes probed by the near limb as a 
function of time. The dashed lines indicate ±1 a uncertainty in the path. The time axis is reversed to be consistent with (a). The occultation midtime 
and immersion/emersion half-light times are indicated. 


tion shadow relative to the observer (the shadow velocity). 
The longitudes and latitudes probed on Triton by each 
chord at the time of half-light are given in Table V (derived 
from Triton’s pole and prime meridian from Seidelmann 
1992). The astrometric solution is depicted in Fig. 5, 
with a globe of Triton as seen from Earth at the event 
time. 


The perpendicular distance between the extreme occul- 
tation chords spans 2013 km on Triton (74% of the globe; 
from the Lick Trl48B chord 706 km north of Triton’s 
center to the WIRO chord 1307 km south of Triton’s cen- 
ter). The mean separation between chords is ~400 km, 
which is go 3d spatial coverage for a small-body occulta- 
tion. 
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TABLE II 

Geodetic Observatory Coordinates from GPS (WGS84) and Half-Light Times 


Site 

Longitude 

a(Longitude) 

(arcsec) 

Latitude 

(»:•:") 

o(Latitude) 

(arcsec) 

Altitude 

(m) 

°alt. 

(m) 

Half-light times UT (mm:ss) after 
1995 08 14 07 hr 
immersion emersion 

Timing error 

(s) 

IRTF 

-155:28:20. 

2.0 

19:49:34. 

3.0 

4182. 

75. 

33:16.29 

35:09.81 

0.08 

KAO 

(immersion) 

-157:30:50. 

1.0 

32:36:10. 

1.0 

13004. 

30. 

38:37.00 

40:17.96 

0.14 

KAO 

(emersion) 

-157:42:55. 

1.0 

32:32:28. 

1.0 

12949. 

30. 

38:37.00 

40:17.96 

0.14 

Lick 

(Trl48A) 

-121:38:38.9 

0.2 

37:20:17.7 

0.2 

1241. 

7. 

36:41.97 

37:56.42 

0.21 

Lick 

(Trl48B) 

-121:38:38.9 

0.2 

37:20:17.7 

0.2 

1241. 

7. 

31:08.60 

32:46.86 

0.79 

Lowell 

-111:32:10.6 

0.2 

35:05:48.7 

0.3 

2209. 

5. 

30:28.78 

32:21.91 

0.36 

WIRO 

-105:58:38.1 

0.6 

41:05:50.0 

0.6 

2956. 

21. 

36:13.67 

37:01.18 

0.74 


TABLE III 


Astrometric Solution from Lightcurve Timings for Trl48 


Parameter 

Circular figure 

Elliptical figure 

Equatorial radius, km a 

1427.9 ±3.5 

1432.6 ±5.7 a 

Ellipticity 

0 . 

-0.029 ±0.016 

Position angle, deg^ 

0 . 

70.3 ±10.1 

/oA.*™ 

2330 ±3 

23 1 6 ± 9 

?oA- km 

-12 ±6 

23 ±22 

/oB> 1™ 

-5157±6 

-5 146 ±9 

goB.km 

-3385 ± 29 

-3431 ±38 

Trl48A-Trl48B separation, arcsec 

0.3869 ± 0.0012 

0.3874 ±0.0011 

Position angle of Trl48B from Trl48A, deg c 

65.7510.16 

65. 16 ±0.14 

Sum of squared residuals 

1504 

707 

Degrees of freedom (dof) 

7 

5 

Weighted RMS residual per dof 

14.7 

11.9 


a This is the equatorial radius at half light in the shadow plane which is one scale height smaller 
than the equatorial radius in the planet plane due to refraction. For the elliptical figure, the 
semiminor axis is given ( — 19.7 ± 10.1 degrees from the equatorial radius). 

b position angle of semimajor axis measured from /, positive in the direction of g to /, which 
corresponds to North through East on the sky. 

c position angle of Trl48B relative to Trl48 measured from North, positive in the direction of 
North to East. 
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TABLE IV 

Adopted Astrometric Solution for Trl48 


Site 

Star 

Impact parameter 
km 

Immersion 
residual, km 

Emersion 
residual, km 

Velocity 

km/s 

I RTF 

Trl48B 

144.5 ±27.3 

0.8 

-0.1 

25.03 

KAO 

Trl48A 

680.2 ± 5.4 

-2.5 

5.0 

24.80 

Lick 

Trl48A 

1086.6 ±5.4 

0.8 

-3.8 

24.95 

Lick 

Trl48B 

702.8 ±27.3 

-13.8 

20.7 

24.95 

Lowell 

Trl48B 

261.0 ±27.3 

-2.1 

-13.8 

24.92 

WIRO 

Trl48A 

1304.1 ±5.4 

-5.9 

-24.7 

24.88 


VII. THE LIGHTCURVE MODEL 

Elliot and Young (1992) describe the processes that af- 
fect the observed flux of an occulted star, and we expand 
their development to explicitly include the far-limb flux 
contribution (which was not necessary for Pluto). For more 
discussion of the stellar occultation technique as applied 
to planetary atmospheres, see the review article by Elliot 
and Olkin (1996). The basic process is depicted in Fig. 6. 
Parallel rays of starlight are incident on the planet from 
the left, with the observer located in the shadow plane a 
distance D from the planet plane. Both the shadow plane 
and planet plane are perpendicular to the incident light 
rays (which are parallel to the .v axis). The coordinate r is 
the radius in the planet plane measured from the center 
of the occulting body. The coordinate y is the observer’s 
location in the shadow plane; the absolute value of y is 
the shadow-plane radius p. As seen in Fig. 6, the observer 
can receive flux from two regions of the limb: the near and 
the far. The near-limb flux received by the observer at y 


has not crossed the x axis, whereas the far-limb flux has 
been refracted across the x axis. 

Starlight with a closest approach to the planet r is bent 
by the refraction angle 6{r ), measured negative toward the 
planet. At any point in the atmosphere, the light ray is a 
distance r' from the center of the planet where r' 2 = r 2 + 
x~. We have assumed the refraction angle is small (6 <§ 1). 

The planet-plane radius r is related to the shadow-plane 
radius p by the refraction angle 6 and the distance from 
the observer to the planet D: 

y(r) = r + D9(r), 

P( r ) = \y{r)\- ^ 

The starlight at the observer is dimmed by (1) spreading 
of the light due to refraction and (2) extinction (scattering 
and/or absorption) of the starlight with line-of-sight optical 
depth T 0 b s . The observed flux is also affected by geometric 
focusing, which is most significant near the center of the 


TABLE V 

Triton Longitude and Latitude Probed at Half-Light by Trl48 


Site 

Star 

Immersion 

Longitude (deg.) Latitude (deg.) 

Emersion 
Longitude (deg.) 

Latitude (deg.) 

KAO 

Tr60 

74.0 

-22.5 

231.1 

3.9 

IRTF 

Trl48B 

81.9 

-6.0 

270.9 

13.5 

KAO 

Trl48A 

110.4 

-27.0 

244.6 

-8.9 

Lick 

TrI48A 

132.4 

-36.5 

227.5 

-22.0 

Lick 

Trl48B 

63.6 

9.8 

291.7 

27.7 

Lowell 

Trl48B 

78.2 

-2.8 

274.5 

16.3 

WIRO 

Trl48A 

153.0 

-40.6 

212.8 

-30.5 
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FIG. 5. Triton’s globe with the six occultation chords. The three 
northernmost chords were occuitations of the fainter star. Trl48B. and 
the three southernmost chords were occuitations of Trl48A. The outer 
circle describes the fitted half-light radius (1427.9 km) in the shadow 
plane, about one scale height smaller than the half-light radius in the 
planet plane due to refraction (Baum and Code 1953). The scale bar 
indicates 500 km or 24 milliarcsec. 


shadow; the starlight is focused by the ratio of the planet- 
plane radius to the shadow-plane radius rlp(r) as the light 
at r is compacted into a smaller circumference (radius p) 
due to refraction. From Elliot and Young (1992) we have 
the following equation for the flux received by the observer 
from the near limb of a spherical planet £(r): 



dr 

dp{r) 


e~ T obs< r \ 


(3) 


Using a radiative-thermal conduction model of atmo- 
spheric structure for Triton (Strobel et al. 1996; hereinafter 
S96), we find the product DO is greater than the surface 
radius for the deepest levels probed by the stellar occuita- 
tions. Therefore, the observer will record flux from both 
the near and far limbs unless there is extinction sufficient 
to completely absorb and/or scatter the flux. 

Now that we have expressions for £(r) and y(r), we can 
formulate the observed flux from both limbs. The observed 
flux from the near limb as a function of y is £(y) - £(r[ v]). 
For a spherically symmetric atmosphere, the far-limb flux 
at y is equal to the near-limb flux at —y, and the observed 
flux from both limbs </>(p) is 


If the body is oblate and p is small enough, the observer 
is within the evolute of a central flash (Elliot et al. 1977) 
and there are more contributions to the observed flux (flux 
is received from all four perpendicular limb points). We 
do not address the flux within the evolute, since none of 
the chords probed within the evolute. 

The fitted parameters in the lightcurve model include 
the full-flux level — the combined light from the occulted 
star and occulting body. For the Trl48A occuitations, this 
also includes the light from Trl48B, since they are unre- 
solved in ground-based imaging. Similarly, for the Trl48B 
events, the full-flux level includes the light from Trl48A. 
Another model parameter, the background-flux level, cor- 
responds to the full-flax level without the contribution 
from the occulted star. A linear slope centered on the 
midtime of the occultation is included in the lightcurve 
model. These parameters are defined in Elliot and Young 
(1992); the units of the background and full levels are sec' 1 
and for the slope it is sec' 2 . 

The remaining parameters are the half-light radius, r h , 
and energy ratio, A h (gravitational potential energy :kT). 
The ratio of these two parameters (r h /A h ) equals the scale 
height at the half-light radius. //(r h ), for an isothermal 
atmosphere. To calculate the model lightcurve, we need 



FIG. 6. Stellar occultation by a planetary atmosphere. Starlight inci- 
dent from the left encounters a planetary atmosphere where refraction 
bends the light rays. The refracted light, dimmed by the spreading of the 
rays, is observed in the shadow plane. Light that has crossed the x axis 
(from where the observer is) constitutes the far-limb flux. The resulting 
occultation lightcurve is seen on the left. The observed signal is the sum 
of all flux at the observer’s location. At some locations this comes from 
one ray; closer to the center of the shadow path there will be two rays 
(from the near and far limbs): and even closer the observer is within the 
evolute of the central flash and there will be more rays contributing 
to the observed flux. (Reproduced, with permission, from Elliot and 
Olkin, 1996). 


4>{p) = C(y) + f(-y). 


(4) 
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TABLE VI 

Parameters for Conversion of Fit Results 


Parameter 

Value 

Reference 

Geocentric distance to Triton (1993/07/10), AU 

29.169 

DE21 1, NEP016 

Geocentric distance to Triton (1995/08/14), AU 

29.269 

DH21 1, NEP016 

Triton mass, kg 

2.1398 10 22 

Anderson et al. 1992 

Molecular weight of N 2 , amu 

28.01 


N 2 refractivity, v$tp at 0.7 Jim 

2.98x10““* 

Peck and Khanna 1966 

N 2 refractivity, vgjp at 2.2 

2.94x10““* 

Peck and Khanna 1966 


to know the distance from the observer to Triton which 
we have approximated by the geocentric distance (see Ta- 
ble VI). We assume the atmosphere is composed entirely 
of N 2 . A small fraction of CH 4 and CO will not significantly 
change the refractivity at standard temperature and pres- 
sure. To convert from lightcurve model parameters (r h , 
A h ) to atmospheric temperature and pressure, we used the 
mass of Triton and the mean molecular weight of the atmo- 
sphere given in Table VI. The refractivity was corrected 
for the actual temperature and pressure. 


VIII. LIGHTCURVE FIT RESULTS 

Examination of the lightcurves revealed that the two 
IRTF lightcurves have a central peak offset from the mid- 
time of the occultation. Furthermore, one can see the slope 
leading up to the central flash is different from immersion 
to emersion (Fig. 7), indicating that the atmosphere of 


Triton near where the central flash is probing is not homog- 
enous. The isothermal model fits assume that the atmo- 
sphere is symmetric, and this assumption is most critical 
near the midtime of the occultation when a wide range of 
Triton’s limb is being probed. Thus, the asymmetry in 
the central-flash lightcurves motivated us to perform two 
different fits to the data. Initially the whole lightcurve was 
fit; then we repeated the fits to a subsection of the data 
that had all points within 20 % of the midtime (where 100 % 
is measured between half-light times) removed. The partial 
lightcurve fits reduce the impact of any asymmetry probed 
as the starlight swept around Triton’s limb. The fitted pa- 
rameters and their formal errors from these fits are given 
in Table VII. The upper row of each parameter corre- 
sponds to the full lightcurve fit, whereas the lower number 
comes from the partial lightcurve fits, which excluded data 
nearest the midtime of the event because this section of 
the lightcurve is most sensitive to asymmetries in the atmo- 
spheric structure as it sweeps over a large section of Tri- 



FIG. 7. IRTF (vis) lightcurve with the midtime from the best-fitting isothermal m idel. The lightcurve demonstrates asymmetry between 
immersion and emersion in two ways: (1) the peak of the central flash is offset from the midtime determined by the midpoint of the fitted half- 
light times, and (2) the slope of the central flash is steeper on emersion than immersion. Bee. .use of this asvmmetry , our adopted fits of the occultation 
lightcurves do not include data within 20% of the occultation midtime. 
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Equivalent-Isothermal Model Parameters from Fits of the Lightcurves' 
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For the Tr60 light curve the midtime is in seconds past 1993/07/10 22:30:00 UTC; for Trl48 it is seconds past 
Fixed (not fitted) parameters. 

Number of degrees of freedom in fit. 
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ton’s limb. Thus, only the main drop and recovery were 
fit in the partial lightcurve fits to minimize the impact of 
our assumption of a spherically symmetric atmosphere on 
the fit results. 

Also given in Table VII is the minimum radius probed 
in the fitted data (r min ). The radius levels probed by the 
individual lightcurves depend on the impact parameter 
of the chord and the atmospheric model describing the 
refractivity of Triton's atmosphere. Chords probing closer 
to the shadow center will probe the atmosphere to a lower 
altitude. At the midtime (for the isothermal atmospheric 
model in Table VII), the IRTF (vis) lightcurve probed 
a minimum radius of 1372 km (~20-km altitude). The 
minimum radius level probed by the IRTF infrared 
lightcurve differs from the visible dataset by ~6 km due 
to the different isothermal models derived from fitting the 
two lightcurves. 

The occultation lightcurves were normalized with the 
fitted background level, slope, and full-flux levels and are 
displayed in Fig. 8. For some of the datasets, only a section 
of the fitted data is displayed. 

From the fit results, we can determine the intensity ratio 
of Trl48A to Trl48B using the Lick Observatory occulta- 
tion lightcurves which recorded both events. From the dif- 
ference in the fitted full-flux level and the background 
level, the intensity ratio of the two stars is 4.19 ± 0.22 
from the unfiltered SNAPSHOT system (similar to an R 
magnitude, Dunham et al. 1991) and 4.44 ± 0.46 from 
the LIRC2 observations which used a K' filter. From the 
magnitude of the blended star V = 12.96 ± 0.03 (personal 
communication from B. Argyle at the CAMC) and K = 
10.98 ± 0.11 (Olkin 1996). we get the individual magnitudes 
of Trl48A ( V = 13.19 ± 0.03, V - K = 1.99 ± 0.12) and 
Trl48B ( V = 14.73 ± 0.05, V - K = 1.93 ± 0.15) if we 
assume the intensity ratio in R equals the intensity ratio 
in V, and that the ratio in K' is the same as that in K. The 
V-K colors are consistent with a spectral type of K0 (Allen 
1973), but if these stars were a physical double, we would 
expect Trl48B to have a later spectral type than Trl48A 
because it is fainter. 

Next, we address extinction using the fit results. We 
determined the zero-flux level of the occultation lightcurve 
and compared it to the fitted value using the intensity ratio 
from the SNAPSHOT system and post-event photometry 
recorded just after the occultation with the same instru- 
ment (to establish the intensity ratio of the blended star 
to Triton). The two methods of photometric calibration 
agreed within their errors and to within 1% of the back- 
ground level (Olkin 1996), which is consistent with there 
being no significant extinction in the region of Triton’s 
atmosphere probed by Lick (since the isothermal model 
assumes a clear atmosphere ). From Voyager data, an inter- 
mittent haze is seen on Triton up to an altitude of 30 
km and condensation clouds exist below 10-km altitude 


(Pollack et al. 1990). Since most of the chords did not probe 
much below 30 km and none probed below 10 km, the 
lack of discernible extinction is consistent with Voyager ob- 
servations. 

For the infrared Trl48B lightcurve recorded at the 1.1- 
m telescope at Lick, the half-light radius and energy ratio 
had to be fixed (not fitted parameters) because of the low 
SNR. They were fixed at the weighted average of the half- 
light radius and energy ratio from the (eight) other Trl48A 
and Trl48B occultation lightcurves. For these data (Lick 
IR lightcurve of Trl48B), only the signal levels were fit, 
and then they were used to determine the intensity ratio 
of the two stars (described above). For the WIRO 
lightcurve, the background level and full-flux level were 
fixed and not fitted parameters. Fits of the WIRO data 
with the photometric parameters free gave inconsistent 
and unrealistic results, which is not surprising due to the 
non-photometric conditions during the observations. The 
full-flux level was determined from post-event photometry. 
To determine the background level, we used the intensity 
ratio of the Trl48B to Trl48A described above (4 19 ± 
0 . 22 ). 

We also performed separate immersion and emersion 
fits to the highest SNR lightcurves (see Table VIII), again 
excluding data within 20% of the midtime. With these 
separate fits, we can look at variations in atmospheric pa- 
rameters with latitude. Figure 9 shows that the derived 
equivalent isothermal temperature, pressure, and half- 
light radius do not vary significantly with latitude. 

Our adopted results for the equivalent-isothermal atmo- 
sphere model of Triton are given in Table IX. They are 
the partial lightcurve fit results of Table VII, with the error 
in the adopted impact parameter factored into the errors 
in the fitted parameters and derived quantities. 

IX. DISCUSSION 

In this section, we address the issues raised in the Intro- 
duction: (1 \ Now that we have probed the altitude range 
20—100 km of Triton’s atmosphere, how do the results 
compare w th published models? (2) Can we detect atmo- 
spheric pre ;sure change with varying insolation? (3) Does 
the structui e of Triton’s lower atmosphere varv with loca- 
tion? We s art with the question of spatial variation. We 
expect Triton s atmosphere to be spherically symmetric 
because the vapor-pressure equilibrium should isother- 
malize the surface by sublimation and condensation, and 
at the higher altitudes (575 km) the atmospheric densities 
at ingress and egress from the UVS occultation data are 
consistent 1 3 within 5% which is consistent with a spheri- 
cally symmetric atmosphere, but not conclusive because 
of the smal sample set (two points on the limb). 

At the altitudes probed by the main drop and recovery 
of the lightcurves, (-30-100 km), the atmosphere shows 
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FIG. 8. Normalized Triton lightcurves and best-fitting model. The partial lightcurve fits did not include data within 20 % of the occultation 
midtime, as indicated by the bar. 


no significant asymmetry (Tables VIII, IX, Fig. 9) in scale 
height and half-light radius. Therefore, the atmospheric 
structure at 1400 km (where the atmospheric pressure is 
best determined) and at 1448.1 ± 1.8 km (the mean half- 
light radius) does not change with location; however, at 
the lowest altitude probed by these occultations (—20 km), 
the IRTF central-flash lightcurves are not consistent with 


a symmetric atmosphere model. More on the inhomogene- 
ity of the atmosphere at this level will be presented in 
Elliot et al (1997). For the rest of the discussion, we use 
the higher altitude levels where we have established that 
the atmosphere is spherically symmetric to within the lim- 
ited sampling and measurement error. 

Next, we address temporal variability in Triton’s atmo- 
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FIG. 8 — Continued 


sphere. Atmospheric change has been predicted to occur 
on Triton because of seasonal variations in insolation caus- 
ing the surface temperature to change (Spencer 1990, Han- 
sen and Paige 1992, Spencer and Moore 1992). The temper- 
ature of the surface ices should become isothermal by 
sublimation and conduction. A change in ice temperature 
will cause a change in atmospheric pressure due to vapor- 
pressure equilibrium. At the time of these observations, it 
was midsummer in Triton’s southern hemisphere as the 


subsolar point on Triton moved in latitude from -48° in 
1993 to -49° in 1995. Unfortunately, due to the single- 
chord natuie of the observations, the uncertainty in the 
pressure determined from the 1993 occultation is too large 
to provide a useful constraint on the pressure change be- 
tween 1993 and 1995. 

We cannot directly compare these occultation results to 
either the v oyager 2 RSS or UVS results because they 
probe a cliff* trent altitude range; however, we can compare 
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Seconds Past 1995/08/14 07:30:00 UT 

FIG. 8 — Continued 


the results of the lightcurve fits to radiative-thermal con- 
duction models based on Voyager observations, such as 
the S96 model. The lightcurve fits give an equivalent-iso- 
thermal description of Triton’s atmospheric structure, but 
the S96 model does not. Therefore, to compare the two we 
need to convert the radiative-thermal conduction model to 
an equivalent isothermal one. This was done by con- 
structing synthetic occultation lightcurves from the S96 


model and then fitting the two-limb isothermal model to 
the synthetic lightcurves (see Table X). Three different 
synthetic lightcurves were constructed from the S96 model 
to simulate the three highest SNR lightcurves (using the 
impact parameter, shadow plane velocity, and integration 
time of the data sets). The atmospheric pressure at 1400 
km derived from the equivalent-isothermal S96 model is 
almost a factor of 2 less than that from the occultation 
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FIG. 8 — Continued 


data. The factor of 2 difference in atmospheric pressure 
derived from the occultation data and the equivalent-iso- 
thermal S96 pressure can be explained three different ways: 
(1) the atmospheric pressure has increased between 1989 
and 1995, (2) the surface radius is larger than the value 
adopted in the isothermal model, and (3) the atmospheric 
structure at the microbar level is not properly modeled. 
Any combination of these could explain the observational 
results, each of which will be explained in turn. Most pub- 
lished seasonal change models do not predict an increase 
in atmospheric pressure of a factor of two or more from 
1989 to 1995, but the basic parameters (e.g., thermal inertia 


and N 2 inventory) that determine the seasonal transport 
of volatiles on Triton are not well known. All the Spencer 
and Moore (1992) models that show an increase in atmo- 
spheric pressure for 1989-1995 have a high thermal inertia. 

The next explanation we consider for the pressure differ- 
ence between the model atmosphere and the occultation 
data is an underestimate of the surface radius. To reconcile 
the factor of 1.8 difference in pressure at 1400 km and a 
scale height of —20 km, the surface radius would have to 
be increased from 1352 to 1363 km. The surface radius of 
Triton was determined by Voyager 2 several different ways 
(with different data sets and different methods). From limb 


TABLE VIII 

Model Parameters from Fits of the Lightcurves (Immersion and Emersion Only) 


Parameter 

Light Curve 

IRTF (vis) 
Trl48B 

KAO 

TrI48A 

Lick (vis) 
Trl48A 

Pmin (1™) 


144.5 ±27.3 

680.215.4 

1086.6 ±5.4 

rh(km) 

immersion 

1446 . 1 !^ 

L46.513.7 

1450.215.1 


emersion 

1449. 

1447.614.1 

1448.115.2 


immersion 

73.614.1 

<9.116.0 

70.317.7 


emersion 

75.714.3 

<9.217.0 

65.717.0 

H(rh) (km) 

immersion 

19.611.1 

:0.9± 1.8 

20.612.2 


emersion 

19.111.1 

<0.9 ±2.1 

22.012.3 

T iso (K) 

immersion 

45.212.5 

48.114.2 

47.215.2 


emersion 

43.912.5 

48.014.8 

50.615.4 

p 1400 (Mbar) 

immersion 

l a <+0.26 
L36 -0.21 

1 40^^ 
1 4U -0.25 

1 as+0-46 
L68 -0.40 


emersion 

■ 3 -0.25 

1 48 +0-30 
J 45 -0.29 

1 54 * 0-40 
1 -^-0.35 



TRITON’S ATMOSPHERIC STRUCTURE 


197 


a b 




c 



FIG. 9. Equivalent-isothermal temperature, atmospheric pressure at 1400-km radius, and half-light radius as a function of latitude on Triton. 
The separate immersion (solid circles) and emersion (open circles) fits of the highest SNR lightcurves probe Triton's atmosphere at six locations 
(from Table VIII). There is no significant trend in these quantities with latitude. This indicates that the asymmetry seen in the central flash does 
not extend to the fitted region of the lightcurves. Also, there is no signature of the 3 % ellipticity of Triton’s atmosphere determined in the elliptical- 
figure astrometric solution. This would be a change of ~40 km (or two scale heights) in the half-light radius. The data scatter less than expected 
for the error bars because a correlated error (the uncertainty in the impact parameter) is included in the error estimate. 


TABLE IX 

Adopted Equivalent-Isothermal Atmospheric Parameters 


Parameter 

KAO 

Tr60 

IRTF (vis) 
Trl48B 

IRTF (IR) 
Trl48B 

KAO 

Trl48A 

Lick (vis) 
Trl48A 

Lick (vis) 
Trl48B 

Lowell 

Trl48B 

WIRO 

Trl48A 

Pmin (1™) 

359 ± 133 

144.5 ±27.3 

144.5127.3 

680.21 5.4 

1086.615.4 

702.8127.3 

261.0127.3 

1304 1 5.4 

Hi (km) 

1451.0^89 

57.5!™ 

1447/7 -2 8 
74.613.0 

1442.9;*;® 

65.017.3 

1447.013.3 

69.114.5 

1449.214.6 

67.915.2 

1452.8!j 7 ; 0 

60.3118.2 

1453.5! 79 
74.91 14.8 

1461.118.1 

51.816.5 

H(Hi) (km) 

25.2 ±3.0 

19.410.8 

22.212.5 

20.911.4 

21.411.6 

24.117.2 

19.413.8 

28.213.6 

T IS 0 (K) 

57 7 + 8-2 
3 // -7.6 

44.51 1.8 

51.315.8 

48.113.2 

48.913.7 

54.9116.5 

44.218.7 

63.618.0 

p 1400 ()* bar ) 

1 8 +7 9 
- 1.2 

14 - 0.2 

1.210.3 

1.410.2 

i £+0-4 
■ 6 -0.3 

1 q+1 8 
1 *- 1.2 

2 O * 0 9 
zu ~ 0.8 

o o+0.9 
2 - 8 - 0.8 
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TABLE X 

Comparison of Adopted Isothermal Atmospheric Parameters and 


Atmospheric Models ' 1 

Parameter 

Light Curve 

IRTF (vis) 
Trl48B 

KAO 

Trl48A 

Lick (vis) 
Trl48A 

Pmin (l® 1 ) 


144.5 ±27.3 

680.2 ± 5.4 

1086.6 ±5.4 

Hi (km) 

data 

1447.7^2 g 

1447.0 ±3.3 

1449.2 ±4.6 


S96 

1434.3^2', 

1434.3 ±2.6 

1434.4 ±4.1 


data 

74.6 ±3.0 

691 ±4.5 

67.9 ±5.2 


S96 

66 - 7 -0.2 

66 7 ±0.3 

66.0 ±0.4 

H(r h ) (km) 

data 

19.4 ±0.8 

20 9 ±1.4 

21.4 ± 1.6 


S96 

21.49 ±0.05 

21.51 ±0.05 

21.72 ±0.07 

Tjso 0Q 

data 

44.4 ±1.8 

48.1 ±3.2 

48.9 ±3.7 


S96 

50.3 ± 0.3 

50.3 ±0.3 

50.8 ± 0.4 

1400 (Mbar) 

data 

1 4 403 
1 4 -0.2 

1.4 ±0.2 

\ 4 
'6-0.3 


S96 

0.8 ±0.1 

0.1 ±0.1 

a 5 -o.i 


a Upper number from adopted fit; lower from fit of S96 model. The error in the parameters in the 
S96 model fits arise mainly from the uncertainty in the impact parameter. 


measurements of Voyager imaging data, the surface radius 
was determined to be 1350 ± 5 km (Smith et al 1989). 
This was refined by more extensive analysis of the imaging 
data by Davies et al (1991), who established a surface 
radius of 1352.6 ± 2.4 km, Gurrola et al (1992) reported 
that the chord length of the RSS occultation gives a radius 
of 1359.6 km with an uncertainty of 0.6A, where A is the 
cross-track uncertainty of the spacecraft. Coupling this 
with the UVS solar occultation gives 1351 ± 6 km, while 
coupling this with the UVS stellar occultation gives 
1342 ± 8 km. Clearly, there is some uncertainty in the 
surface radius of Triton, but nevertheless most solutions 
would not be consistent with a radius of 1363 km. An 
underestimate of Triton’s surface radius may not be the 
sole explanation for the pressure difference, although it 
could be part of the answer. 

The last explanation for the larger surface pressure is 
that the radiative-thermal conduction model does not cor- 
rectly describe Triton’s atmosphere in the microbar region. 
The S96 model is based on RSS measurements that probed 
near the surface and UVS solar occultation data that 
probed much higher altitudes (500 km). The structure be- 
tween these two regions is dominated by heat conduction 
down from the ionosphere and radiative transfer by CO 
and CH 4 . If this model does not accurately describe the 
region probed by the occultation, we would expect the 
occultation data and radiative-conductive models to dis- 
agree (as they do). 

We have other evidence that the physics in the region 


probed by the occultation may not be correct: the strong 
thermal gradient predicted by a radiative-thermal conduc- 
tion model (Strobel et al 1996) may be present in our 
highest SNR data, but shifted to a colder temperature. The 
highest SNR data (IRTF vis) were fit using the power-law 
thermal gradient model adapted from Elliot and Young 
(1992) to include the far-limb flux contribution and resulted 
in a thermal gradient of 0.29 ± 0.25 K/km, (see Fig. 10) 
again the data within 20 % of the midtime were excluded 
from the fit. This gradient is consistent with that from the 
S96 model at 1400 km (0,14 K/km) and the Krasnopolsky 
et al (1993) model (~0.2 K/km). The temperature at 1400 
km from t ie thermal-gradient model fit is 41.5 ± 5.2 K — 
colder than the S96 model (52.4 K) and consistent with 
the Krasnopolsky et al (1993) model. 

Consider the isothermal fit results to the occultation data 
in a different way. From the adopted fits (Table IX), the 
weighted average pressure at 1400-km radius in 1995 is 
1.4 ± 0.1 /^bar with an equivalent-isothermal temperature 
of 47 ± 1 K. If we assume the surface pressure is still 14 
/xbar (Gurrola 1995) as measured in 1989, then using the 
atmospheric pressure at 1400-km radius of 1.4 ^bar, we 
would expect an equivalent-isothermal temperature of 53 
K. This temperature is higher than either the measured 
temperature of the surface ice 38 K (Tryka et al 1993) or 
the equivalent isothermal temperature derived from fitting 
the highest SNR occultation lightcurves (47 ± 1 K). We 
can reconcile these two temperatures (53 K and 47 K) if 
one of our assumptions is incorrect. The two assumptions 
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0 10 20 30 40 50 60 70 

Temperature (K) 

FIG. 10. Atmospheric temperature from the occultation data, two 
atmospheric models, and the Voyager RSS experiment. The RSS radius 
is from Tyler etal (1989) and the RSS equivalent-isothermal temperature 
is from Gurrola (1995). The equivalent isothermal temperatures as deter- 
mined from model fits are indicated by filled symbols. The isothermal 
temperatures refer to the range of atmosphere probed by the occultation 
and are indicated where the minimum fractional pressure and thermal 
gradient errors are near 1400 km. The ler temperature profile from the 
power-law thermal gradient model fit of the IRTF (vis) data is indicated 
by the bars. 


were (1) the surface pressure is still 14 jubar, and (2) the 
surface radius is 1352 km. We can reduce the 53 K by 
either increasing the atmospheric pressure or the surface 
radius. We have already addressed these possibilities. 

X. CONCLUSIONS 

We have fit an isothermal model to the 10 occultation 
lightcurves that constitute the complete set of ground- 
based stellar occultation data to date. Comparisons across 
the separate immersion and emersion lightcurve fits show 
no significant trend of temperature, pressure, or half-light 
radius with latitude. The consistency of the fitted half-light 
radius (1448.0 ±1.8 km) and of the atmospheric pressure 
at 1400 km (1.4 ± 0.1 fibar) indicates that the atmosphere 
is spherical at those levels within the errors and sampling 
limitations; however, we must keep in mind that the cen- 
tral-flash observations indicate the atmosphere is not 
spherically symmetric at the deepest levels probed, —20- 
km altitude. 

We have also compared these occultation results to radi- 


ative-thermal conduction models of Triton’s atmosphere 
and have shown that the atmospheric pressure and temper- 
ature gradient predicted by radiative-thermal conduction 
models do not match that derived from the occultation 
data. The pressure at 1400 km is about a factor of 2 larger 
than that in the S96 model at the same radius, and the 
temperature derived from the occultation is consistently 
colder. The difference in atmospheric pressure could be 
due to seasonal change, but given the difference between 
the model and observed temperature, we believe at least 
some of the difference is a result of an inadequate model 
in the altitude range 30-100 km. 

We look forward to regularly probing Triton’s atmo- 
sphere with future stellar occultations to investigate the 
effects of seasonal change on Triton’s tenuous atmosphere. 
Multi-chord occultation observations will allow us to better 
constrain seasonal change in atmospheric pressure. On 
1997 July 18, Triton is predicted to occult the star Trl76 
(McDonald and Elliot 1995). The star is brighter than any 
of the previously observed Triton occultation stars (V = 
12.53 ± 0.02, J = 11.2 ± 0.1, K = 10.5 ± 0.2; the V 
magnitude is courtesy of R. Stone at the USNO). It is 
unfortunate that the KAO has been shut down because 
the predicted track of this occultation is mostly over the 
Pacific Ocean. 
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as Tycho 651672 and GSC632 1-0 1030) was successfully observed 
with the Hubble Space Telescope (HST) in daylight over the 
northwest Pacific Ocean; Astrometer 3 of the Fine Guidance Sensors 
(FGS) was used to record the event 5 . Details of the HST data are 
■ given in Table 1 along with information about other observing 
stations; the immersion and emersion data (disappearance and 
reappearance of the star) are shown in Fig. 1. A central flash (the 
focusing of light rays by Triton’s atmosphere 6,7 ) was recorded, but 
will be presented and analysed elsewhere. 

We modelled the HST light curve with a standard small-planet 
model that allows for a power-law temperature gradient* (Table 2). 
The background from dark counts and Triton (determined by the 
FGS in September and adjusted for Tritons different distance) was 
subtracted, and the remainder was divided by the flux from the star 
(also determined by the FGS in September) so that the full range of 
stellar flux corresponded to values between 0.0 and 1.0. In the light- 
curve model fits, the zero level was fixed, but the full-scale signal 
from the star was a free parameter. The difference of the fitted values 
from 1.0 shows that our calibration error is only a few tenths of one 
per cent. 

From fitting the entire light curve, the closest-approach distance 
between the centre of Triton’s shadow and the HST was determined 
to be 224 ± 4 km (first column of results in Table 2). This value 
places the shadow somewhat further north than predicted, but it is 
consistent with no detectable occultation at our Oahu station 
(Table 1). Using the closest-approach distance determined from 
the entire light curve, we fitted the main immersion and emersion 
sections of the light curve both separately and together (next three 
columns of results in Table 2). As a test of the self-consistency of our 
light-curve models, we fixed the atmospheric model parameters 
(“half-light radius”, “lambda at half-light,” and the “thermal- 
gradient exponent”; see ref. 8) at their values determined from 
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Triton, Neptune’s largest moon, has been predicted to undergo 
significant seasonal changes that would reveal themselves as 
changes in its mean frost temperature 1 ' 5 . But whether this 
temperature should at the present time be increasing, decreasing 
or constant depends on a number of parameters (such as the 
thermal properties of the surface, and frost migration patterns) 
that are unknown. Here we report observations of a recent stellar 
occultation by Triton which, when combined with earlier results, 
show that Triton has undergone a period of global warming since 
1989. Our most conservative estimates of the rate of temperature 
and surface-pressure increase during this period imply that the 
atmosphere is doubling in bulk every 10 years — significantly 
faster than predicted by any published frost model for Triton 2,3 . 
Our result suggests that permanent polar caps on Triton play a 
dominant role in regulating seasonal atmospheric changes. Simi- 
lar processes should also be active on Pluto. 

The 4 November 1997 occultation of the star Trl80 (also known 



Figure 1 Triton occultation light curves from the HST Data before and after the 
occultation were used to establish the modulation of the signal due to the 
astrometer scan, which was then removed from the entire data set The zero and 
full-scale stellar flux levels were established with photometric data from an earlier 
FGS visit to these objects on 11 September 1997 At that time, Triton's magnitude 
as observed by the FGS was 13.4 and the magnitude of TrlSO was 10.6. The FGS 
data have been averaged at 1.0 s; a, immersion data (filled circles} and light-curve 
model fit (line); b, emersion data (filled circles) and light-curve model fit (line; see 
fits in Table 2). The zero point of the abscissa is arbitrary. The light-curve model 8 
used a power-law thermal gradient; residuals from the fit in a are shown in c, and 
the residuals from the fit in b are shown in d. The r.m.s. residual along the full 
signal is -0.0022 for a 1-s average, and is the result of photon noise. The 
remaining residuals that occur when the star is partially occulted are the result 
of unmodelled structure in Triton's atmosphere and rarely exceed 0.01. The effect 
of these on our determination of the 1,400-km pressure can be estimated by 
differences among the fits in Table 1. The omitted central portion of the light curve 
(see text) corresponds to radii <1,400 km (altitudes <49 km); the atmosphere 
below that does not affect the 1,400-km pressure. 
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the combined immersion and emersion fit (second to last column in 
Table 2), and then fit the entire light curve again (last column in 
Table 2). 

Although the agreement is not perfect between all model fits for 
the atmospheric parameters at a radius of 1,400 km, it is very good. 
In particular, the close agreement between the immersion and 
emersion atmospheric parameters for two widely separated loca- 
tions on Triton supports the idea that the sublimation and con- 
densation of nitrogen maintains the surface frost at the same 
temperature. The pressure at 1,400 km, 2.3 ± 0.1 p,bar, is derived 
from the fit to the immersion and emersion data (the error bar 
includes the systematic differences between the fits in Table 1). This 
is significantly greater than the 1,400-km pressure of 1.4 ± 0.1 fxbar 
measured with stellar occultations in 1995 9 and the value of 
0.8 ±0.1 fLbar extrapolated from a surface pressure of 
14 ± lpibar measured by Voyager in 1989 10,11 . The pressure at 
1,400 km has certainly increased between 1995 and 1997, based 
on the two sets of stellar occultation data; the temperature at 
1,400 km has also increased from 47 ± 1 K (ref. 9) to 
50.3 ± 0.5 K. As these temperatures are consistent with the 1,400- 
km temperature predicted by atmospheric models based on Voyager 
data 12,13 , a surface-pressure increase between the 1989 and more 
recent measurements is the most likely explanation for the differ- 
ence (rather than using an inappropriate model for extrapolating 
the pressure from the surface to 1,400 km). 


This surface-pressure increase implies a temperature increase of 
the surface frost, as the principal constituent of the tenuous atmos- 
phere, nitrogen, is presumed to be in vapour- pressure equilibrium 
with the surface frost. The surface pressures and corresponding 
equilibrium temperatures for N 2 are plotted in Fig. 2, where two 
cases are shown. 

In the first case, we have extrapolated the occultation pressures 
from a radius of 1,400km to the surface radius of 1,352 km (an 
average of published values 14,15 ) using a multiplicative factor of 1 7.5, 
which is valid under the assumption that the shape of the thermal 
profile 12 has not significantly changed since the Voyager encounter 
(filled symbols in Fig. 2). We justify this assumption by the 
consistency of the atmospheric temperature at 1,400 km (Table 2) 
with the models 12,13 . These measurements show a steady increase in 
pressure and surface-ffost temperature. At this rate the atmosphere 
has tripled in bulk since the time of the Voyager encounter, which 
would require ^*0.9% of the total solar energy incident on Triton 
during the intervening time to sublime the frost In the second case, 
we make a more conservative estimate of the surface pressure 
derived from the stellar occultation data by assuming an isothermal ! 
atmosphere below 1,400 km (open symbols in Fig. 2). Even for this ! 
case, the atmosphere has nearly doubled in bulk since the time of 
Voyager. J 

Triton’s atmosphere is supported by surface frosts, which at any 
instant are expected to be at the same temperature everywhere, due j 


Table 1 Observation* 


Site 

Telescope aperture 
(m) 

Detector 

Filter 

Recording interval 
(mon 4 Nov. 1997) 

Integration time 
(s) 

Aperture 

(arcsec) 

S/N* 

HSTt 

2.5 

4 PMTs $ 

F583WII 

04:59:32-05:24:44 

0.025 

5x5 

630 ; 

LCC4 


CCD 

Nonell 

05:10:55-05:27:35 

0.125 

15 x 15 

33 1 


*S/N refers to the signal-to-noise ratio of the occultation light curve (when the star is unocculted) over a time interval that corresponds to 20 km (approximately a scale height) of relative 
motion between the telescope and Triton's shadow. 

tHST was located (inJ2000 coordinates, kmrelative to the centre of the Earth) at (x,y\z) =«( -3,906, -5,752, 761 ) for immersion, and (x,y,Z) = ( - 1,896, -6.464. 1.806) for emersion. Centre- 
line observations with NASA's Lear Jet Observatory were attempted, but failed due to target acquisition and tracking problems. Simultaneous optical and infrared observations wrth NASA's i 
Infrared Telescope Facility on Mauna Kea (MKO) and optical observations with a 0.4-m portable telescope on Tern Island were foiled by douds; other groups had similar bad weather on the j 
University of Hawaii 2.2-m telescope on MKO and on Kauai (more details are provided at http://occultmitedu/tr180plans/). | 

t The four photomultiplier tubes (PMTs) of Astrometer 3 of the Fine Guidance Sensors were used in "TRANS" mode. We generated the light curve by adding the signals from all four PMTs 5 . j 
4 The telescope at Leeward Community College (LCC) is located at longitude -157" 59' 6* E. latitude 2T23* 30" N, and an altitude of 10 m. The data are a continuous series of CCD frames ! 
(each 60 by 57 pixels). The light curve was generated from these frames with synthetic aperture photometry of the blended image of Triton and the star. This station recorded data through the j 
time of the event but saw no obvious occultation, which indicates that Triton's shadow path went somewhat north of the prediction (provided at http://occultmit.edu/tri80.html). j 

H The effective wavelength of the F583W filter was 583 nm; the open CCD used at LCC had an effective wavelength of 670 nm. | 


Table 2 Models for the Triton occultation light curve obtained with the HST 




Data selection 



AH 

Immersion* 

Emersion* 

Immer. & emer. 

All 

Model parameters 






Background level 

0.00 

0.00 

0.00 

0.00 

0.00 

Background slope ( 10~ 7 km" ' } 

5.3±0.7 

0.00 

0.00 

2.2 ±0.8 

4.9±0.6 

Star signal 

0.9968 ±0.0003 

0.9957 ±0.0004 

0.9974 ±0.0003 

0.9966 ±0.0002 

0.9971 ±0.0003 

Closest approach timet 

3.828±0.034 

3.808 

3.808 

3.808 ±0.01 5 

361 7 ±0.034 

Half-light radius (km) 

1,456.3±0.6 

1,456.3 ±0.3 

1,456.4±0.2 

1,456.3±0.2 

1,456.3 

Lambda (isothermal) at half-light 

65.3±0.4 

65.5±0.4 

66.4±0.3 

67.4±0.2 

67.4 

Thermal-gradient exponent 

f.9±0.3 

1.7±0.5 

4.4±0.4 

2.9 ±0.3 

2.9 

Minimum centre distance (km) 

224 ±4$ 

224 

224 

224 

219±2 

Derived quantities atr = 1,400 km (altitude 

= 48 km) 





Number density ( 10 u cm^ 3 ) 

3.23±0.09 

3.23±0.06 

3.42 ±0.04 

3.31 ±0.04 

3.31 

Pressure (p.bar) 

2.25±0.06 

2.24±0.04 

2.36±0.03 

2.30±0.03 

2.30 

Temperature (K) 

50.6±0.3 

50.4±0.3 

50.1 ±0.2 

50.3 ±0.2 

50.3 

Temperature gradient (K km -1 ) 

0.07±0.01 

0.06 ±0.02 

0.16 ±0.01 

0.10±0.01 

0.10 

Other derived quantities 






Minimum attitude probed (km)4 

<25 

46 

46 

46 

<26 

Surface pressure (*ibar) 

39±4 

39 ±4 

41 ±4 

40±4 

40 

Fit information 






Degrees of freedom 

1,504 

321 

321 

644 

1,507 

x 2 per degree of freedomll 

17.8 

8.94 

4.22 

6.77 

18.1 


These models are for a pure N 2 atmosphere. Model parameters wrth error bars were fitted: those without were held constant. 

• The immersion data included the interval from 05:01 :03 to 05:06:27 uron 4 Nov. 1997; the emersion interval was 05: 10:06 to 05: 15:30 ut. Immersion probed the atmosphere above -5.1* latitude 
and 10.4° longitude on Triton; emersion probed the atmosphere above -8.4° latitude and 206.4" longitude, 
t Seconds after 05:08:00 ut on 4 Nov. 1997. 

t The minimum centre distance is pnmanly controlled by the amplitude of the central flash, which will be discussed elsewhere. 

4 For the fits to all the data, the ‘minimum altitude probed' refers to the light-curve centre. A surface radius of 1,352 km was used' 4 ’ 6 . 

II x 2 per degree of freedom is much greater than 1.0 for all fits, indicating that the deviations of the atmospheric structure from a power-law thermal gradient 8 (rather than photon noise) are the 
primary source of error (note the behaviour of the residuals in the lower panels of Fig. 1 ). 
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to efficient transfer of latent heat through the atmosphere 16 * 17 . This 
global frost temperature determines atmospheric pressure through 
vapour-pressure equilibrium; thus a frost temperature increase 
from 37.5 K (in 1989) to 39.3 K (in 1997) can be inferred from 
the observed increase in atmospheric pressure (the non-isothermal 
case in Fig. 2). 

The mechanisms that can cause a change in frost temperature 
(and hence surface pressure) are: (1) non-static (migrating) surface- 
frost distribution; (2) changes in the optical properties of the surface 
frost; and (3) changing insolation on a static surface-frost distribu- 
tion. (Although subsurface heat flow — geothermal or seasonal — 
can affect the energy balance, it is unlikely to have varied rapidly in 
the past decade.) 

Simple seasonal models of long-term frost migration, with uni- 
form frost albedo and emissivity 1,3 , show monotonic long-term 
transfer of most of Triton’s nitrogen frost to the poles, which have 
the minimum seasonally averaged insolation. In most of these 
models, the remaining seasonal frost migrates between hemispheres 
in response to the seasonal insolation cycle, and all frosts should be 
almost exhausted from the southern hemisphere during the current 
extreme southern summer. This would result in a rapid decrease in 
atmospheric pressure, the opposite of what our 1995 and 1997 
observations indicate has occurred. 

A substantial decrease of mean frost albedo or emissivity could 
account for the pressure change. For instance, model "F” of Spencer 
and Moore 3 , which assumes unit emissivity, can account for the 
atmospheric change if the mean frost Bond albedo dropped by 0.12 
between 1989 and 1997. Observational 18,19 and theoretical 20,21 work 
suggests that large changes in frost optical properties may indeed 
occur on Triton. 

The frosts total absorbed insolation may have increased, due to 
the increase in the southerly subsolar latitude in the past decade. We 
computed the change in surface pressure that would result solely 
from the change in subsolar latitude from August 1989 (47.8° south) 
to November 1997 (51.0° south), under the assumption that no 
volatile transport had occurred in the intervening 8 years. Five 
different maps of the 1989 N 2 -ice distribution on Triton were 
generated based on the albedos measured by Voyager. The emissiv- 



Flgure 2 Triton's global warming. The surface pressure and the corresponding 
vapour-pressure equilibrium temperature is plotted versus year of measurement. 
The measurement in 1989 refers to the results of an atmospheric mode!’ 2 based 
on several measurements made by Voyager 2. The 1995 data are from stellar 
occupations of Tr148A and Tr148B 9 , and the 1997 data refer to the present work. 
Each occupation datum has been plotted twice along with a corresponding least- 
squares fit: assuming that the atmosphere follows the model of Strobe! et al' 2 
(filled circles, solid line), and making the conservative assumption that the 
atmosphere is isothermal below 1,400 km (open circles, dashed line). For the 
conservative assumption, the pressure increase is 1.4 ± OAiAbaryr"’; it is 
2.8 ± 0.8p.baryr 1 for the other case. In both cases, the pressure is seen to 
steadily increase, indicating that Triton has been undergoing a period of global 
warming. The distance of the points from the lines is larger than their error bars, 
which indicates that the warming of the surface frost may not have occurred 
linearly with time. 


ity of the N 2 ice was taken as a free parameter in each of these frost 
models, and was adjusted such that the temperature of the N 2 was 
equal to 37.5 K (corresponding to 14 pibar surface pressure) in 1989. 
The N 2 -ice thermal balance was then recomputed using the subsolar 
latitude appropriate for Triton in late 1997, under the assumption 
that the distribution of N 2 ice was unchanged (albedos did change 
slightly, due to the illumination-angle dependence of the bolometric 
albedo, and the rotation of some N 2 -ice-covered portions of the 
globe into darkness). Several of these models result in pressure 
increases comparable to the measured increase. Thus, the observed 
global warming of Triton may be due to increased insolation of a 
permanent south polar cap. □ 
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